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Abstract. We present a proposal for a versatile cold-atom-based quantum 
simulator of relativistic fcrmionic theories and topological insulators in arbitrary 
dimensions. The setup consists of a spin-independent optical lattice that traps 
a collection of hyperfine states of the same alkaline atom, to which the different 
degrees of freedom of the field theory to be simulated are then mapped. We 
show that the combination of bi-chromatic optical lattices with Raman transitions 
can allow the engineering of a spin-dependent tunneling of the atoms between 
neighboring lattice sites. These assistcd-hopping processes can be employed 
for the quantum simulation of various interesting models, ranging from non- 
interacting relativistic fermionic theories to topological insulators. We present a 
toolbox for the realization of different types of relativistic lattice fermions, which 
can then be exploited to synthesize the majority of phases in the periodic table 
of topological insulators. 
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1. Introduction 

In a seminal paper published in 1982 [1 , R. P. Feynman discussed in great detail 
the problems connected with the numerical simulation of quantum systems. He 
envisaged a possible solution, the so-called universal quantum simulator, a quantum- 
mechanical version of the usual simulators and computers currently exploited in many 
applications of the "classical" world. If realized, such a device would be able to 
tackle many-body problems with local interactions by using the quantum properties 
of nature itself [5] ■ Interestingly, even without the advent of a fully universal quantum 
computer, the construction of small dedicated devices, also known as purpose-based 
quantum simulators, would already be of significant importance for our understanding 
of quantum physics. The basic idea is to engineer the Hamiltonian of the quantum 
model of interest in a highly-controllable quantum system, and to retrieve all the 
desired information with a measurement of its properties. Many research fields would 
eventually benefit from such devices, such as two- and three-dimensional many-body 
physics, non-equilibrium dynamics, or lattice gauge theories [3]- 

Recently, the scientific community is considering ultra-cold atoms as one of 
the most promising candidates for the realization of a wide variety of dedicated 
quantum simulations [HE]. Indeed, these gases are genuine quantum systems where 
the available experimental techniques offer an impressive degree of control together 
with high-fidelity measurements, thus combining two fundamental requirements for 
a quantum simulator. Among the most recent experimental achievements, we 
would like to mention the observation of Anderson localization in disordered Bose- 
Einstein condensates (BEC) [SllZ], the research on itinerant ferromagnetism with cold 
fcrmions [5], or the reconstruction of the equation of state of fermionic matter in 
extreme conditions, such as in neutron stars [S]. 

An important drawback in the applicability of cold atoms as quantum simulators 
is the difficulty of coupling their spatial degrees of freedom to external magnetic fields. 
This prevents a direct simulation of the quantum Hall physics [10j . the controlled 
observation of whose extraordinary phenomenology would shed new light on quantum 
many-body theory. One way to overcome this problem is to dress the system with 
ingenious laser schemes, which mimic the effect of an external magnetic field, and thus 
allow the neutral atoms to behave as effectively charged particles [IT]. This approach 
led recently to the realization of neutral BECs coupled to external effective magnetic 
and electric fields [121 [13], or even with an effective spin-orbit coupling [14]. More 
generally, the scientific community has now realized that even in presence of an optical 
lattice, dressing cold gases with suitable optical and microwave transitions could 
push the experiments beyond the standard superfluid - Mott insulator transition, and 
significantly widen the spectrum of the models that are currently being simulated |15) . 
The possible applications of such optical-lattice-based quantum simulations are 
numerous and diverse, ranging from the realization of Abelian and non-Abelian static 
gauge fields [TT1 H3 [13 HH1 HH to that of quantum Hall states [HI [2JJ [23 (23 [23] ; horn 
the study of the anomalous quantum Hall effect [551 120] to the quantum spin Hall 
effect [371 HE] ; from three-dimensional topological insulators [2H] , to flat-band physics 
with a non-trivial topological order [30 , or non-Abelian anyons [31]. Recently, a big 
effort has also been put in designing schemes where the exotic effects associated to 
relativistic quasiparticles, such as the Klein tunneling and the Zitterbewegung, arise in 
a controlled table-top experiment [3^l3l?l l3 ^ l3l ^ l3l ^ l3T l l3 ^ l3l?]. 

In this article, we elaborate on the idea of using a spin-independent bi-chromatic 
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Figure 1. Optical superlattice potential of equation ([TJ in the two-dimensional 
case, with parameters Vb = WE r and f = 1. Left: the potential is characterized 
by a square geometry of main minima; in the middle of each link an intermediate 
minimum is also present. Right: if the lattice is deep enough, the spectrum of the 
system features two energy bands whose Wannier functions are localized in the 
main minima and in the secondary minima, as plotted in the figure. 



optical lattice dressed with suitable Raman transitions to simulate interesting non- 
interacting field theories of lattice fermions. We present a concrete proposal to 
create a three-dimensional optical lattice that traps a multi-species atomic gas, and 
to tailor arbitrary spin-dependent hopping operators. We have already shown how 
this setup could break the SU(2) invariance of the hopping rates for spin-1 atoms in 
spin-independent lattices, and how the simulation of systems subjected to three-body 
repulsion could benefit from it [3D]. Here, we extend this idea further, and show that 
the same setup allows for the realization of hopping operators which modify the atomic 
hyperfine state. Combining this trapping scheme with Fermi gases, we show that this 
platform would open a new route towards the simulation of high-energy physics and 
topological insulators. 

This paper is organized as follows: in section [2j we describe qualitatively the idea 
of using an optical superlattice to realize a general hopping operator for a multi-species 
cold gas of alkalis. Further analysis and technical details are given in section [3j where 
we also present some numerical results that support the possibility of controlling a 
spin-flipping tunneling in this platform. The reader not interested in these technical 
details may skip this content without prejudicing the comprehension of the following 
sections. Some final remarks on the proposal are presented in section [4j In section [5j 
we discuss the possible applications of the described scheme, focusing on relativistic 
theories and topological insulators, and trying to give a list of the most interesting 
phenomena which could be explored. Finally, we present our conclusions in section [6] 

2. The Setup and The Idea 

We consider the following atomic three-dimensional optical potential 

V(x) = -V [cos 2 (^-)+^cos 2 (2 9 x J )] , (1) 

j'e{i,2,3} 

where x = (2:1, 2:2, £3), q = 2-k/Xl (Xl is the wavelength of the laser), and where 
Vq, £ > represent the potential amplitudes. The low-energy structure of this potential 
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Figure 2. Sketch of the atomic structure of 40 K: from the electronic structure 
(L is the electronic angular momentum) to the fine structure (J = L + S, where 
S is the electronic spin) to the hyperfine structure (F = J + I, where I is the 
nuclear spin). The latter is drawn in the specific case of an external magnetic 
field present. The last box shows the optical spin-independent potential which 
traps equally all the hyperfine levels. 



is a cubic array of main minima separated by "secondary" minima located in the 
middle of each lattice link (see figure [l}. We note that additional higher-order minima 
are also present, but will not play any role in the phenomena discussed in this article. 
Due to the specific form of the potential in equation ([I]), the Hamiltonian can be 
divided into three independent terms, each one depending on one of the three couples 
of conjugate operators, {2!j,Pi}tei,2 3- Consequently, the Bloch functions of the n- 
th band with energy E n (p), can be written as tp n>p (x) = Ylj (^j)- I* 1 or der to 
discuss the effects occurring on the scale of one lattice site, Wannier functions can be 
introduced for each band 

w„,r(x) = — J e~ lRp V>„, P (x)<ip = w n , Rl (xi) w n ,R 2 (x2)w n ^ 3 (x 3 ). 

Like Bloch functions, Wannier functions belonging to different bands form an 
orthonormal basis, and one can thus expand the Hamiltonian in such a basis. Since 
the Wannier functions are not eigenstates of the Hamiltonian, this expansion leads to a 
Hubbard model describing the tunneling of atoms between neighboring sites, together 
with a local on-site interaction coming from the scattering of the cold gas [41] . 

This setup can be used for the simulation of a lattice field theory, where the field 
operators are identified with the atomic creation-annihilation operators in the Wannier 
basis of the lowest energy band (i.e. the states localized in the main minima of the 
lattice). Conversely, higher energy bands provide auxiliary levels that shall be used 
as a resource to tailor the tunneling processes. The main result of this article is the 
claim that a complicated though not unfeasible combination of current technologies 
leads us to the realization of the following Hamiltonian 

H sys = [Uu]t'tCtt + nC tr' [Mr'rCrr + H.C. (2) 

ru tt' r 

Here, we are considering a multi-species fermionic scenario with many hyperfine levels 
of the same atom: c\ T (c rT ) creates (annihilates) a fermion with hyperfine spin r 
localized in the main minima of the superlattice at r = m\B.i + mi&i + m^a.^, where 
rrij 6 {l...Lj}, Lj stands for the number of lattice sites along the and a., is 

the lattice spacing in the j-th direction. The parameter t v stands for the strength of 
the laser-assisted tunneling in the i> direction, with with v g {ai, a2, a^}, which shall 
be described below. The operators U v describe the tunneling from r to r + i/, and are 
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Figure 3. Sketch of a laser-assisted tunneling induced in the presence of a 
superlattice. Two physical hyperfine states belonging to the F = 9/2 manifold 
are connected via Raman couplings with the intermediate level of an auxiliary 
state belonging to the F = 7/2 manifold. If the coupling is detuned enough, 
the F = 7/2 level can be adiabatically eliminated: no population is left there 
and an effective coupling is engineered between neighboring sites. Left: scheme 
for a spin-preserving (i.e. diagonal) hopping. Right: scheme for a spin-flipping 
hopping. 



a common feature in lattice gauge theories. We have also included an on-site Raman 
term A, of strength f2, that induces a certain transition between the hyperfine states. 
Note that we use Gaussian units and h = 1. The claimed possibility of engineering a 
wide range of hopping operators t/„, together with the state-of the-art control of the 
atomic interaction, makes already our system a versatile quantum simulator of lattice 
field theories. In this manuscript, we focus on non- interacting theories, which can be 
realized either with dilute systems, or by employing Feshbach resonances to lower the 
interaction strength (see e.g. Oil]). We stress that interesting phenomena can also be 
observed in non-interacting gases when additional ingredients are introduced in their 
dynamics, such as disorder or the assisted-hopping processes discussed in this article. 

Let us note that the control of the homogeneous tunneling for a single-species 
atomic gas is straightforward, and would not even require the superlattice (£ = 0) [42 . 
Moving to a many-species case, one runs into the problem that a general hopping 
operator also entails terms flipping the atomic hyperfine spin (simply referred as spin 
in the following), which are not easily engineered. Here, we propose to realize such 
couplings by combining Raman transfers and a bi-chromatic superlattice (£ ^ in 
equation Q). The proposal can be applied to all the alkalis notwithstanding their 
bosonic or fermionic nature. In the following, however, we shall focus in the fermionic 
scenario, which is best explained with the following practical example. 

Let us consider an ultra-cold cloud of non-interacting 40 K atoms in the presence 
of a magnetic field of intensity B. Such field lifts the spin degeneracy within the two 
atomic hyperfine manifolds of the ground state, F — 9/2 and F — 7/2, according to 
the following relations (see also figure [2]): 

- B 9/2,m F = +9F fJ-B B m F E 7 / 2 ,m F = A H f - 9f Hb B m F (3) 

where m F is the projection of the hyperfine spin along the quantization axis defined 
by the magnetic field, /is is the Bohr magneton, g F is the hyperfine Lande Factor, 
and Ahf stands for the hyperfine splitting. These hyperfine levels are all trapped into 
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Table 1. Numerical values of a possible three-dimensional optical bi-chromatic 
superlattice |lj for 40 K used in section [3] for numerical simulations. We 
characterize the properties of the two energy bands, which exhibit Wannier 
functions localized respectively in the main and secondary minima, by listing 
the energy expectation value (Ej) of the most localized Wannier function and the 
bandwidth AEj. Finally, we argue that atoms trapped in optical lattices show 
a hierarchy of typical energies which can be actively exploited for engineering 
non-trivial hopping operators. 
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the same spin-independent optical potential Q. Depending on the lattice theory we 
want to simulate, we select a subset of these hyperfine levels described theoretically 
by creation-annihilation operators in the lattice sites. We then identify such fields 
with the components of the lattice field theory to be simulated. This leads us to 
divide the hyperfine levels into two subsets: the subset of "physically meaningful" 
states, which belong to the hyperfine manifold F = 9/2, and the usually larger subset 
of auxiliary levels that shall be used to assist the tunneling and create the desired 
hopping operator. 

Regarding the hopping operator in equation ([2]), we address each of its matrix 
elements [£/„] T ' T separately. Given a matrix element (i.e. once we have identified 
the initial and final hyperfine levels to be connected by the assisted tunneling), we 
choose an auxiliary level belonging to the hyperfine manifold F = 7/2 trapped in 
the middle of the link. These levels provide intermediate "bus" states that shall be 
used as a resource to assist the tunneling as follows. The couplings between the 
atoms in the main sites, Ri, and the "bus" states, R2, are realized via optical two- 
photon Raman processes transferring a net momentum q t . They have a mathematical 
expression proportional to the overlap integral of the initial and final Wannier 
functions: J ui* 2 Rs (x)e ,q '' x w Iii Ri (x)dx. This integral is not zero because of the term 
e iqt " x , which is of course relevant only if 27r/|q t | is of the order of the lattice spacing. 
Since this regime cannot be achieved with microwave transitions, one is motivated 
to employ two-photon Raman transitions. Interestingly enough, it is possible to 
eliminate adiabatically the intermediate level and obtain an effective four-photon 
coupling between neighboring sites (see figure [3]) . We stress that different matrix 
elements can be engineered at the same time thanks to the magnetic-field splitting of 
the hyperfine levels ([3]): the involved atomic transitions become non-degenerate and 
can be individually addressed with different lasers. Furthermore, the use of coherent 
laser light for the Raman transitions entails the additional advantage of being able 
to deal with complex phases, and thus to realize complex gauge structures at will. 
The realization of the non-diagonal matrix elements requires the lattice to be slightly 
staggered, a technique discussed also in reference [17] . Summarizing, this proposal 
tries to exploit a hierarchy of energies characterizing atomic gases in optical lattices 
in order to assist the tunneling between neighboring sites with controlled adiabatic 
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eliminations (see table [lj. 

The on-site spin-flipping A in equation ([2]) can be performed with standard 
technology based on microwave transitions, or Raman transitions carrying negligible 
momentum. Furthermore, these terms can also be exploited to correct spurious on- 
site couplings which may be induced by the laser scheme. Unfortunately, we note that 
there is no selection rule relying on the polarization properties of the light and the 
hyperfine moment of the atoms, which can be used to realize the different on-site and 
nearest-neighbors spin-flipping processes. The superimposed magnetic field cannot be 
aligned at the same time with the propagation vector of all the three lasers, aligned 
along the three Cartesian axes, which would be the case in which circularly polarized 
light could be exploited to induce controlled transitions. Conversely, we shall rely 
on the different Zeeman-shifted energies to selectively address the different couplings 
between the internal states. 

3. Realization of Spin-Dependent Hopping Operators 

In this technical section, we theoretically and numerically confirm the qualitative 
scheme presented above. We study two simple but important cases: the realization of 
diagonal and non-diagonal hopping operators for a two-species atomic gas. These can 
be considered as the main building blocks needed to realize any tunneling operator 
even in situations with more than two atomic species. 

3.1. Coupling Between Different Hyperfine Manifolds 

The most fundamental ingredient of this proposal is the possibility of using Raman 
processes to induce controlled atomic transitions between different hyperfine states of 
the electronic ground state L = [L is the total electronic angular momentum). These 
transitions are realized with two lasers via adiabatic elimination of the electronically 
excited manifold L = 1. In the following, we address the atomic levels as \L,a,k), 
with a labeling the hyperfine degrees of freedom (see also figure [2] for some insights 
on the internal structure of 40 K), and k the quantum numbers of the center-of-mass 
wavefunction (in our case, the Wannier functions of the optical potential). As discussed 
in [ID], the induced Raman coupling between the state |0, a,k) and |0, a',k') can be 
written as follows: 

&a'k';ak{t) — Sk'k ^a'a 6 (4) 

This expression clearly factorizes the following contributions: 

• the time-dependence of the effective coupling and its effective frequency, which is 
the difference between the frequencies of the two lasers lj — u>i — lj 2 ; 

• the dependence on the center-of-mass degrees of freedom, Sk'k = 
(fc'|e~ l( - P2 ~ Pl " ) x |/c), where pi and P2 are the momenta of the two lasers; 

• the dependence on the initial and final internal states and on the polarization 
properties of light, D a / a , which is a function of the dipole matrix elements between 
the initial (final) state and the excited levels. 

Next, we specify Q to the superlattice setup of section[2j i.e. we will consider Raman 
transitions in presence of lattices characterized by a Wannier function trapped in the 
middle of each link. 
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\k = l) |J6 = 3) 

Figure 4. The "6-level model" used to model the spin-preserving (diagonal) 
hopping of F = 9/2, m F = 9/2. The auxiliary state F = 7/2, m F = 7/2 has 
been chosen. The center-of-mass quantum number is k. Energies are not in scale; 
the orders of magnitude of the parameters are the following: d ~ 10 -4- 100 kHz, 
8 ~ 100 -f 300 kHz and A ~ 1 -4- 10 GHz. We propose to adiabatically eliminate 
the upper manifold and to study the dynamics of the lowest one with an effective 
Hamiltonian H per t BSD ■ 



3.2. Developing an Effective "6-Level Model" 

Let us address the simulation a theory characterized by two-component fields. 
Following the discussion of section [2j we take two states of the F = 9/2 manifold 
of 40 K, for instance 1 9/2; mp = 7/2) and 1 9/2; mp = 9/2), and map them into the 
theory to be simulated. Here and in the following subsections, we discuss the laser- 
assisted hopping in the diagonal case (mp preserved while hopping) and non-diagonal 
case (mp flipped while hopping). 

For the diagonal case, we develop the "6-level model" depicted in figure [4] We 
consider one physically meaningful state, say \F = 9/2, mp = 9/2), and one auxiliary 
state, say \F = 7/2, mp = 7/2). Moreover, we consider different Wannier states for 
each of them, two localized in main sites (k = 1 and 3) and one in the intermediate 
link (k = 2). The model includes the effects of undesired couplings and additional 
levels, and its limitations, together with the approximations on which it relies, will 
be discussed at the end of the paragraph. We can identify the states with the short 
notation \F,k) rather than with the longer previous one \0ak). Below, we give an 
analytical estimate of the population transfer rate, whereas in the next subsections 
we present the numerical time-evolution for physically interesting cases. 

The model is parametrized by six relevant couplings between the different Wannier 
functions Sk>k (see figure [4]), whose properties are listed below. We exploit the 
existence of theorems which assure the possibility, in our case, of considering three real 
and exponentially localized Wannier functions Wj(x),j g {1,2,3} [43]. We write the 
parameters Sk'k factorizing out the space dependence of the coupling e iqt ' x J , where 
Xj is the position of the point around which the Wannier function Wj(x) is localized, 

S k ' k = e**'** J <,(x-x fe , +x fc )e lqt - x ^(x)dx; (5) 

S\,l = / 5*2,2; 5*1,3, 53,1 ~ 0. (6) 

The parameters 5i,i and 52,2 describe two on-site couplings, whereas 5i,2 is the 
coupling between a main site and an intermediately trapped state (see figure |4| . The 
last relation states that couplings between neighboring main sites are negligible. The 
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relation between the other four overlap factors depends on the particular experimental 
situation. In this case, we are interested in the simplest scenario where a single Raman 
transition induces all these couplings, which leads us to 



Si, 



S 2 i = e 2ic lt- x i e i qf( x 2- x i) > 5'* _ e 2»qfXi e iqt-(x3-xi) <jf 



J 2,3- 



In order to make this scheme simpler, we assume q t = 2qx, and thus e jq *'( X2 ~ Xl ) = 
1. As we will argue below, transferring a momentum which does not fulfill this 
requirement is not a problem since the resulting phase can be gauged away. The 
phase 2q t ■ xi can also be put to zero for the moment, since its role only becomes 
important when one needs to give a phase to different matrix elements. In the 
following, we will also consider situations where the coupling between the lattice 
sites 2 and 3 could be induced by lasers propagating in the opposite direction, where 
Si t 2 = 5*2,1 = 5*3,2 = 52,3- Taking these considerations into account, the Hamiltonian 
reads as follows (see figure [4] for the definitions of S, A and u): 

H = d |9/2, 2) (9/2, 2| + (A + d)|7/2, 2) (7/2, 2| + 

+ A(|7/2,l)(7/2,f| + |7/2,3)(7/2,3|) + 



+ Qe- wt [5'i,2(|7/2,2)<9/2,l|- 
+ 5^ (17/2,2X9/2,31 
+ 5 M (|7/2,1)(9/2,1| 
+ 5 2 , 2 |7/2,2)(9/2,2| [ 
Once we apply the unitary transformation 



17/2,1X9/2,21) - 
-|7/2,3X9/2,2|)- 
h 1 7/2,3X9/2,3|) 
+ H.c. 



(7) 



r(t) = exp[id(|9/2,2)(9/2,2| + \7/2, 2)(7/2, 2|) t], 



the three levels |9/2, k) become degenerate. In case the three inequalities \Si jfl\/(6 — 
rf) < 1 are fulfilled, it is possible to use second-order perturbation theory in order to 
develop an effective Hamiltonian describing the dynamics within the sub-manifold we 
are interested in, namely 



Hp er tJQ 



!5i,i| 2 , |5 lj2 | 2 



S-d 5 
[|9/2,l)(9/2,l| + |9/2,3)(9/2,3[ 
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, |5i, 2 | 2 
'S -2d 



o2 
°1,2 
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1 



5 
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|9/2,2)(9/2,2| 
- H.c. 

1 > idt 



5 -2d 



'2,2 5l,2 



id I 



S-d 

1 1 

2 \S^d + S 

■ [|9/2,2)(9/2,l| + |9/2,2)(9/2,3|] + H.c. (8) 

Remarkably enough, this Hamiltonian leads to the desired transfer rate of population 
from level |9/2,1) to |9/2,3), and viceversa. The main contribution is the direct 
coupling 



J 



13 e 



|5i, 2 | 2 fi 2 



,i20. 



arg 5i 



(9) 
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A second contribution, which in our system will prove to be not-negligible, comes from 
a sort of "adiabatic elimination" of the level |9/2, 2), namely 

_ .(2) = (9/2, 3|ff pert |9/2, 2) (9/2, 2|g pert |9/2, 1) 

13 (9/2,2|J? pert |9/2,2) - (9/2, l|# pert |9/2, 1) + cf 1 j 

Accordingly, we have derived the desired effective Hamiltonian where the Raman lasers 
assist the hopping of the physically meaningful F = 9/2 levels, after the auxiliary 
F = 7/2 bus states have been adiabatically eliminated. In the following sections, we 
shall address the range of validity of the approximations leading to this Hamiltonian, 
and compare it with the exact numerical investigation of the initial Hamiltonian Q. 

We want to stress here that even if the integrals in the definition ([5| of the Skk' 
can be complex numbers, this does not have any physical influence on this proposal. 
Indeed, even if the effective coupling between neighboring main sites — J was complex, 
its spatially uniform phase can be gauged away with a space-dependent unitary 
transformation (even in the case of periodic boundary conditions). Conversely, the 
non-uniform phase coming from the e jqt ' Xfc factor, which arises when q t is not parallel 
to the direction of the hopping it assists, cannot be gauged away even in presence of 
open boundary conditions. Such a phase, which is not related to the fact that the 
integrals in ([5| are complex, can be used to simulate an external uniform magnetic 
field j!6l I17j . Finally, we underline that in our setup, where the tunneling along each 
axis is induced by lasers propagating parallel to the axis itself, both complex phases 
can be gauged away. In order to simulate a magnetic field, therefore, one should move 
slightly away from this configuration and engineer a Raman coupling whose effective 
transmitted momentum does not run parallel to the links of the lattice. We will not 
consider this situation in this article because the models of interest in section [5] do not 
require such space-dependent phase. 



3.3. Range of validity of the "6-Level Model" 

The presented "6-level model" strongly relies on two approximations: 

(i) considering the bands of the lattice as being flat; 

(ii) neglecting delocalized higher-energy free states. 

If these approximations are not justified for a given experimental configuration, 
spurious population transfers to next-neighboring sites would arise. 

The approximation ^ is required to fulfill the core idea of the proposal, namely 
the adiabatic elimination of the intermediate level. This is demonstrated with a model 
which considers only a subset of the Hilbert space spanned by the real eigenstates of the 
Hamiltonian (Bloch functions), considering just three of their linear combinations (the 
Wannier functions u?fe=i(x), u>2(x) and W3 (x)). This is equivalent to approximating 
the dispersion laws of the band as being flat, neglecting thus possible curvature effects, 
and is legitimated as long as the width of the band is much smaller than the detuning of 
the transition 6 — d. In case the degeneracy of the Bloch functions cannot be assumed, 
all the Bloch functions should be considered in order to quantitatively estimate the 
spurious effects cited above. In general, this issue sets a trade-off for the relative depth 
£ of the secondary lattice in ([!]): on one hand, a shallow lattice (£ < 1) is desirable 
because the Wannier function of the intermediate minimum Wk=2 (x) is not strongly 
localised and laser- induced transitions are favored (|5i,2| ~ | <Si,i |)- On the other 
hand, the more the wavefunction is delocalized, the more the band bends, eventually 
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Figure 5. Left: sketch of the scheme proposed for the realization of a diagonal 
hopping operator (Energies are not in scale). Raman coupling 1 (2) connects 
the \F = 9/2, m F = 9/2} (\F = 9/2, m F = 7/2}) states to their auxiliary state. 
Detuning allows independent control of the hopping rates. Right: exact time- 
evolution of the "6-levels model" Q, showing the coherent population transfer 
between sites 1 and 3 of the spin state \F = 9/2, mp = 9/2}. The parameters 
used are listed in table [5] The inset shows the maximal populations of the six 
considered levels labelled \F,k) as in |7| and shows that only a small fraction of 
the population is lost in auxiliary levels. 



becoming parabolic at k = with a bandwidth comparable to the detuning. In our 
numerical simulations we consider £ = 1, which is a reasonable middle- way. 

Regarding the issue higher-energy bands could become important in the 
presence of intense Raman transitions f2 and large detunings d — d, which couple 
them to the lowest-band states. The presented analytical and numerical studies do 
not take into account these effects since they consider only three Wannier functions 
and effectively only two bands, even though including bands with localized Wannier 
functions would just imply a renormalization of the numerical coefficients Sk,k'- A 
different problem is the case of high-energy strongly parabolic bands, whose Wannier 
functions are not strongly localized. The effect of such states is not considered by 
our model, which is that of spreading population among many next- and further- 
neighboring main sites. From an experimental point of view, we expect a trade-off to 
arise between a large detuning regime, allowing powerful lasers and strong effective 
couplings with noisy spurious population transfers, and a small detuning one, with 
clean but small couplings. 

3.4. Diagonal Hopping Operator 

We now explicitly study the possibility of realising a diagonal tunneling operator. We 
numerically simulate the Hamiltonian Q with a simple Runge-Kutta algorithm. We 
did not include in the simulation hypcrfine states different from \F = 9/2, mp = 9/2) 
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Table 2. Parameters used for the numerical simulation of the diagonal hopping 
in subsection |3.4| We list the numerical values of all the main parameters 
characterizing the atomic transitions and the Raman couplings. The first Raman 
coupling induces the hopping of the F = 9/2, mp = 9/2 whereas the second one 
addresses the F = 9/2, mp = 7/2 (such states were however not considered in 
the simulation). 



Level: \F, mp, k) 


Energy 


Parameters 


|9/2; 9/2; 1) 


g F ^ B Bm F 


A HF 


1.285 GHz 


|9/2;9/2;2) 


g F fj, B Bm F + d 


fi F B 


40 MHz 


|9/2; 9/2; 3> 


g F fj, B Bm F 


d 


69.160 kHz 


|7/2; 7/2; 1) 


A H f - g F ^i B Bm F 


5*1,1 


0.46 


|7/2; 7/2; 2) 


Ahf — 9F\J-BBm F + d 


5*1,2 


0.07 + i0.13 


|7/2; 7/2; 3) 


A H f - g F ii B Bm F 


5*2,2 


0.16 



# Raman 


Q 




1 


49.5 kHz 


^|7/2;7/2;2) 


— ^|9/2;9/2;l) 


- 300 kHz 


2 


49.5 kHz 


^|7/2;5/2;2) 


— -E|9/2;7/2;2) 


- 300 kHz 



J 



'13 J if | Estimated T 
176 Hz 17 Hz I 0.018 s 



Numerical T 
0.017 s 



Table 3. Parameters used for the numerical simulation of the non-diagonal 
hopping in subsection |3.5| We list the numerical values of the main parameters 
characterizing atomic transitions and Raman couplings. The reported frequencies 
of the Raman couplings are approximate because some additional tuning is needed 
to compensate the different Stark shift for states with k = 1 and k = 3 arising 
due to Raman dressing in presence of staggering. Perfect matching the atomic 
transition becomes difficult and imperfections are responsible for the not clean 
population transfer in figure [6] Larger staggering values would help. 



Level: 


F, m F , k) 


Energy 


Parameters 


|9/2; 9/2; 1), 


|9/2;7/2; 1) 


g F ^ B Bm F 


A HF 


1.285 GHz 


|9/2;9/2;2), 


|9/2;7/2;2) 


g F [i B Bm F + d 




40 MHz 


|9/2; 9/2;3>, 
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g F [i B Bm F + 15 kHz 


d 
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# Raman 



49.5 kHz 
49.5 kHz 
49.5 kHz 
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|9/2;7/2;l) 
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300 kHz 
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Figure 6. Left: sketch of the scheme proposed for the realization of a non- 
diagonal hopping operator (Energies are not in scale). Raman couplings 1, 2 
and 3 connect the \F = 9/2, m F = 9/2} and \F = 9/2, m F = 7/2} states to the 
auxiliary states. Detuning allows independent control of the hopping rates. Right: 
exact time-evolution of the "12-levels model" introduced in subsection 13.51 and 
generalizing jTJ. We show the coherent population transfer between sites 1 and 3 
of the spin state \F = 9/2, m F = 9/2} and \F = 9/2, m F = 7/2}. The parameters 
used are listed in table [2] The inset shows the maximal populations of the 
twelve considered levels labelled \F,m F ,k) and shows that only a fraction of 
the population is lost in auxiliary levels. 



and \F = 7 /2,mp — 7 /2) because they are strongly detuned from those we are 
considering. However, for completeness, we include the presence of a second Raman 
coupling which would be needed to induce the hopping of \F = 9/2; mp = 7/2) and 
check that it is unimportant. 

We show in figure [5] the numerical results. The realistic parameters used in this 
simulation are listed in table [2] The population is coherently transferred between two 
neighboring levels and only a negligible fraction is lost in auxiliary states. Regarding 
the validity of the "6-levels model" , for the lattice considered here, the bandwidths 
of the two bands are respectively 0.2 kHz and 9.1 kHz, which should be compared 
with the considered detuning of 300 kHz. In these and the following simulations, the 
employed numerical values have only an illustrative purpose and other regimes could 
be considered. 

Therefore, these results confirm the plausibility of our scheme to induce a laser- 
assisted tunneling between the atoms sitting in the main minima of the optical lattice. 
To make the simulation toolbox reacher, we now address the possibility to control a 
spin-dependent hopping process. 
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3.5. Non-Diagonal Hopping Operator 

In order to study the realisation of the non-diagonal hopping operator, we consider 
an enlarged 12-levels model, which is a generalization of the previous one taking into 
account more hyperfine states. We want now to transfer population between the 
manifolds \F — 9/2, mjr = 9/2) and \F — 9/2, TOf — 7/2) and consider as auxiliary 
states \F = 7/2, m F = 7/2) and \F = 7/2 1 m F = 5/2). 

A big issue which must be solved to engineer such a hopping is the arousal of 
undesired spin-flipping terms induced by the laser. In this paper we consider the 
possibility of staggering the lattice with an additional optical field, in order to lift 
the degeneracy between the different sites of the optical lattice, in the same fashion 
of |17j . Such a staggering can be done also in three dimensions since the cubic lattice 
is bipartite, and we consider staggering values of 10 — 15 kHz. 

Figure [6] sketches the experimental scheme we have in mind and shows the exact 
time evolution of the population transfer between the two levels of F = 9/2 in two 
neighboring sites. Interestingly enough, we show a flip of the Zeeman spin during the 
tunneling process, and thus obtain the promised spin-dependent hopping operator. 
The parameters of the simulation can be found in table [3] 

4. From a Spin-Dependent Hopping Operator to a Quantum Simulator 

In the previous section, we discussed how the superlattice geometry could be used to 
create non-trivial hopping operators on each link. Here we want to assemble these 
ingredients and discuss how to use them to engineer interesting quantum simulator in 
arbitrary dimensions, considering the advantages and disadvantages of the proposal. 

First of all, we stress that the lasers needed to engineer the hopping along one 
direction must transfer momentum along that same direction. Therefore, just by 
controlling the beam propagation directions, we can tailor different tunneling operators 
along each axis. This is an important feature which will be largely exploited in the 
proposals listed in section[5] This kind of directionality selection rule is also responsible 
for avoiding the population of higher-order minima which do not lie on the edges of the 
unit lattice cell. Since they are not connected to the main minima by a line parallel 
to a cartesian axis, we do not consider any momentum transfer along such direction, 
and thus the formal orthogonality of the Wannier functions localized in those minima 
is never lifted. 

Due to the very general formulation of the superlattice potential, the setup is 
well-suited also to work in two and one dimensions. Moreover in ID it is possible to 
align the magnetic field splitting the hyperfine sublevels with the optical lattice. In 
this case, it is possible to use specific light polarization to selectively couple different 
atomic levels. As a short-term goal, it would be very interesting to understand what 
is the most interesting physics which could be simulated in a one-dimensional system, 
where the presence of more symmetries could lower the experimental intricacies. 
We partially address this question in section [5j where we argue that several one- 
dimensional topological phases could be realized. In a three- (two-)dimensional case, 
with the magnetic field aligned with the diagonal of the cube (square), the polarization 
of the light cannot be exploited, only energy-based selection rules are reliable. 

Unfortunately, energy-based selection rules are not enough to prevent, in the case 
of spin-flipping operators, spurious on-site spin-flipping couplings. We proposed to 
solve this issue by staggering the optical lattice, i.e. lifting the degeneracy of the 
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lattice sites of some tens of kHz (see also the discussion in reference [17]). One should 
also mention that the diagonal hopping operators can still be engineered in presence 
of such staggering, with the only additional issue of using two Raman couplings (as 
in the non-diagonal case) to match the energy difference between sites. 

Each of the matrix elements of the hopping operators is realized via an effective 
four-photon process. This means that the spin of the atom can be flipped of at most 
| Ami? | = 4: a careful analysis is needed in case one is interested in simulating a theory 
with more than 4 fields, because some hopping matrix elements might become avoided. 

Finally, the description given in this proposal is essentially at the single-particle 
level where no many-body effects have been used. As a consequence, the proposal 
works both for bosons and fermions, which is a valuable result. The proposals of 
section [5] assume the possibility of switching off atomic interactions with Feshbach 
resonances. Nonetheless, it would be of the utmost interest to study how interactions 
could fit into this picture, and to explore the rich variety of models that could emerge. 
We leave this problem for future work. 

We now conclude the part of the article devoted to the description of the 
experimental setup. We believe to have provided the relevant results supporting the 
initial claim that it was indeed possible to realize a system whose low-energy structure 
is described by Hamiltonian 

5. Applications of the Quantum Simulator 

In this section, we would like to discuss the use of the described setup as a quantum 
simulator. We will address a range of lattice field theories for relativistic fermions [44 , 
and explore exotic phases of matter known as topological insulators [35] ■ In general, 
the task of a purpose-based quantum simulator is to realize a system described by an 
effective Hamiltonian H e s that reproduces faithfully the properties of the model to be 
simulated. In our case, this model corresponds to relativistic lattice fermions -ffrei, or 
topological insulators H top . The resource to accomplish such a task is the microscopic 
control over the supcrlattice setup, which we have argued previously to be described 
by the Hamiltonian rewritten here for reading convenience, 

H sys = ^ X! tljC l+uT' [Uu]t'tCtt + ^ Qc\ T , [A] r / T C rr + H.c. 

rv rr' r 

The main objective in the following subsections is to control and manipulate: 

(i) the optical lattice dimension D; 

(ii) the tunneling strengths t v ; 

(iii) the spin-dependent hopping operators U v ; 

(iv) the on-site Raman transitions Q, A; 

such that the Hamiltonian of equation ^ simulates the desired physics, namely 

H sys ({t u , U TV , A, ft}) — > H c ff PS -ffrel, H top . 

From a condensed-matter perspective, exotic phases are frequently associated 
to strongly-correlated regimes and many-body interactions. There are however 
distinguished exceptions to this paradigm, such as graphene [46] and topological 
insulators |45| , where quadratic fermionic Hamiltonians contain a wealth of non-trivial 
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phenomena. In the case of graphene, a two-dimensional layer of graphite, the low- 
energy carriers can be described by emerging relativistic fermions without mass. On 
the other hand, topological insulators are exotic holographic phases with an insulating 
bulk, and a peculiar boundary that hosts robust conducting modes protected by 
topology arguments. In both cases, the transport properties of the material differ 
significantly from the standard solid-state theory. In the subsections below, we show 
how the Hamiltonian ^ serves as a versatile simulator of these two interesting phases 
of matter. 



5.1. The Zoo of Relativistic Lattice Fermions 

The properties of a relativistic spin- 1/2 fermion with mass m are described by the 
famous Dirac Hamiltonian |47) 

H = J dr *(r) t i/ DI *(r), H m = cot ■ p + mc 2 /3; (11) 

where a u ,{3 are the so-called Dirac matrices fulfilling a Clifford algebra, {a^,^} = 
2<5 i//i , {a v , (3} — 0, and c stands for the speed of light. Here, ^(r) is the TV^-component 
fcrmionic field operator, where Np = 2 for one and two spatial dimensions, and 
Nd = 4 for three spatial dimensions. Our objective now is to construct an effective 
Hamiltonian starting from equation ^ that closely resembles the relativistic field 



theory in equation (11). The underlying setup consists of a gas of ultracold 40 K 
atoms, which is a non-relativistic system; nonetheless we can design it as a quantum 
simulator of relativistic particles by exploiting the quantum statistics and a peculiar 
engineerable Fermi surface. 



Naive Massless or Massive Dirac Fermions. The idea is to engineer translationally 
invariant hopping operators, U v = e^ 1 " 4 " , according to the SU(A r £) ) group, where N D 
has been defined above. For the particular choices specified in table |4j one finds that 
the Hamiltonian in equation ^ in momentum space becomes 

H = ^2 *k I 2t " cos $ v cos kvl + 2tlJ sin sin ku0lv * k ' ( 12 ) 

kGBZ \ v ) 

where is a multicomponent fermi operator that contains the different No hyperfinc 
levels involved in the simulation, and k is defined within the first Brillouin zone 
BZ. One readily observes that there are certain regimes, the so-called 7r-flux phases 
(j> v = 7r/2, where the energy spectrum develops Nd = 2 D degeneracy points K d where 
the energy bands touch e(K<j) = 0. Around these points K<j = (d x Tr, d y iT, d z ir), where 
d v E {0, 1} is a binary variable, the low-energy excitations of the 40 K Fermi gas are 
described by the effective Hamiltonian 

^ff^EE^fP^^i^tPd), J ffg I ( Pd )-ca d - Pd , (13) 

d Pd 

where p d = k — K d represents the momentum around the degeneracy points, 
(a d )„ = (— l) d "a v are the Dirac matrices listed in table |4j and c = 2t x = 2t y — 2t z 
is the Fermi velocity that plays the role of an effective speed of light. Therefore, the 
Fermi surface of the half-filled gas consists of a set of isolated points, the so-called Dirac 
points, and the low-energy excitations around those points behave according to the 



Hamiltonian of massless Dirac fermions in equation ( 13 1. As occurs with graphene [46], 
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Table 4. Quantum simulator of naive Dirac fermions. Each translationally- 
invariant hopping operator U v is expressed in terms of Pauli matrices {cr x , a y , cr z }, 
and a set of dimensionless fluxes {01,02,03}- We also list a particular 
representation of the Dirac matrices a u ,fi for different spatial dimensions d, 
together with the important symmetries of the Hamiltonian that depend on T. 



D 




u &2 


u* 3 


1 


gi0io-m 






2 




e i<t>2<Ty 




3 






gi03(T z (g)cr z 



D 



0~7. ® o„ 



<J X ®h 



-iot x a v ot z 



or other cold-atom systems [33J [33J HS1 13H [3H [37] , one finds a table-top experiment 
of relativistic quantum field theories. 

Let us note that we obtain an even number of relativistic-fermion species, each 
located around a different Dirac point (i.e. M\ = 2 for one dimension, A/2 = 4 
for two dimensions, and A3 = 8 for three dimensions). This doubling of fermionic 
species is a well-known phenomena in lattice gauge theories |44j . where the fermions 
in equation (13 1 would correspond to the so-called naive Dirac fermions [48j . As 
predicted by the Nielsen-Ninomiya theorem [33], this doubling cannot be avoided 
without breaking an underlying symmetry: 

• for D odd {H^ 1} Ti} = is an involution known as chiral symmetry; 

• for D even r\ [-Hdi( — Pd)l T2 = -ffpi(Pd) ^ s an antiunitary symmetry known as 
time-reversal (see table |4|. 

According to these results, we have a versatile quantum simulator of massless 
Dirac fermions in any spatial dimension D = 1, 2, 3. In D — 1, 2 they coincide with 
the famous Weyl fermions, whereas in D = 3 they contain a couple of Weyl fermions 
with opposed helicities. Note that this scheme can also be extended so as to simulate 
exotic Weyl fermions of any arbitrary spin s |50) . Besides, our quantum simulator also 
allows us to make these fermions massive, thus reaching the desired Hamiltonian in 
equation pTj ). The idea is to control the on-site Raman transitions such that A = /? 
listed in table[4j In such case, the Rabi frequency plays the role of the mass mc 2 = 2f2, 
and the effective Hamiltonian in equation ( 13 1 becomes 

ff&(p d ) =ca d -p d +mc 2 (3. (14) 

Therefore, this quantum simulator can explore both the non-relativistic and the ultra- 
relativistic limits of the theory. 



Wilson and Kaplan Fermions. From a lattice gauge theory perspective, the 
additional fermions around Kj 7^ are spurious doublers that modify the physics at 
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long wavelengths. A partial solution is to give the doublers a very large mass mK d c 2 , 
so that they effectively decouple from the low-energy physics of the Dirac fermion 
at Kd = 0, namely TOK d ^> tok - We must find a way to engineer a momentum- 
dependent mass that differs from the global on-site Raman mass discussed above. 
By combining the laser-assisted tunneling listed in table |1J with the additional terms 



U„ — ie lv "' 3 [29], the momentum-space Hamiltonian H in equation (12) becomes 
H + H, where 



H = "J^ I 2i v cos ip v sin k v I — 2t v sin cp v cos k„(3 

keBZ \ v / 



*k, (15) 



where t v are the additional laser-assisted tunneling strengths. Once more, for the 
7r-flux phases tp v = ir/2, the effective Hamiltonian in equation (|11| is modified into 



^m(Pd) = ca d • p d + m Kd c 2 /3, m Kd = m - ^(-l) d "m^, (16) 

where to„c 2 = 2t v . Let us emphasize that our quantum simulator provides a complete 
control over the different masses, since m depends on the on-site Raman transition 
strengths, whereas m u depends on the assisted-hopping strength, and thus on the 
laser power. In particular, when these parameters fulfill to„ — m (i.e. m x = m for 
D = 1, m x + m y = to for D = 2, and m x + m y + m 2 = to for Z) = 3), we obtain a single 
massless Dirac fermion at Kd = 0, whereas the remaining doublers have been boosted 
to much higher energies. We thus achieve a quantum simulator of the so-called Wilson 
fermions of any spatial dimension D — 1, 2, 3 |51) . 

We note that this decoupling between a single massless Dirac fermion and its 
doublers is not in conflict with the Nielsen-Ninomiya theorem since the introduced 
mass terms explicitly break the aforementioned symmetries. This is particularly 
important in odd dimensions, where the theory does not preserve chiral symmetry, 
a fundamental concept in the standard model classifying right/left-handed particles 
Ti^J = To preserve such symmetries, the concept of Kaplan fermions arises |52) . 
namely massless Dirac fermions bound to a lower-dimensional domain wall located 
at where the Wilson mass gets inverted TOK d = —\m\ + 2|m|#(rj_ — rj_). Since 
we have a complete experimental access to the parameters of the Wilson mass, it is 
also possible to tune them such that to„ > m, and thus the mass mK d < gets 
inverted, and one gets a lower-dimensional massless fermion bound to the region where 
this mass inversion takes place. 

Let us close this subsection by underlying the versatility of our setup as a quantum 
simulator of a diverse set of relativistic lattice fermions. Not only can we implement 
massless Dirac fermions of any dimensionality, thus exploring their connection to Weyl 
fermions, but we can also control their mass. This leads us to the concept of massive 
Dirac fermions, and the notorious Wilson and Kaplan fermions dealing with the 
fermion doubling problem. Interestingly enough, the physics behind these high-energy 
particles is intimately related to the materials known as topological insulators [45 , 
which are the subject of the following subsection. In fact, it is always possible to find 
a Kaplan- fermion representative within each class of topological insulators [S3] . 

5.2. A Toolbox for Topological Insulators 

Topological insulators correspond to fermionic gapped phases of matter that are 
insulating in the bulk but allow a robust transport along the boundaries [45]. This 
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Table 5. Periodic table of topological insulators. The underlying quadratic 
Hamiltonians H = ^Z a /3 /3 + H.c, where a, /3 represent the lattice 

sites and the internal states of the fermion, can be classified according to 
the fundamental symmetries of time-reversal T, charge conjugation C, and the 
combination of both S = TC. The values T = 0, C = 0, S = 0, are assigned 
to Hamiltonians that break the symmetry, whereas T = ±1, C = ±1, 5 = 1 
correspond to symmetry-preserving Hamiltonians, where T 2 = ±1, C 2 = ±1, 
cS 2 = +1. There are six possible combinations for non-interacting fermionic 
Hamiltonians (and another four for pairing fermionic Hamiltonians), which lead 
to the classes listed in the first column. For each dimension d, there are three 
possible topological insulators among all these classes, and they are classified 
according to the integer Z or binary Z2 nature of a certain topological invariant. 
In the column labelled QS, we list the particular instances that can be simulated 
with our super-lattice based quantum simulator. 
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Unitary 











AIII 


Chiral unitary 
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yes 
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2Z 
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yes 



robustness is due to the existence of gapless edge excitations which are protected 
against disorder by topological arguments (i.e. they avoid Anderson localization [54], 
and thus transport charge even in the presence of strong disorder) . The paradigmatic 
example of a topological insulator is the integer quantum Hall effect (IQHE) [IS], a 
two-dimensional electron gas subjected to a strong magnetic held that displays a robust 
quantization of the transverse conductivity a xy — n ■ e 2 / h, where n G Z is related to 
the topological invariant known as the Chern number [56] . In this case, chiral electrons 
bound to the one-dimensional edges of the sample avoid back-scattering processes, and 
are thus immune to disorder [57] . Remarkably enough, the IQHE is only one instance 
of a large list of topological insulating phases. Each class can be characterized by a 
set of discrete fundamental symmetries |58| . and a certain topological invariant, see 
table [5] Note that we have excluded the topological superfluids from this table, since 
their quantum simulation would require a pairing mechanism, and thus goes beyond 
the scope of this work |591 160) . We now discuss how our quantum simulator can 
reproduce the properties of many of these fascinating phases of matter following two 
possible strategies. 
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Bottom-Up Approach. In this case, one designs the ultracold-atom Hamiltonian so 
that it simulates a particular model belonging to the desired class of topological 
insulators. Therefore, a different experiment would be required for each class-oriented 
simulator. Two representative, yet reasonably simple examples are the Su-Schrieffer- 
Hegger model of polyacetilene [ST], which is related to the D = 1 BDI topological 
insulator, or the tt— flux phase of the fermionic Creutz ladder [62], which is related 
to the D — 1 AIII topological insulator. The former can be simulated by using a 
one-component Fermi gas in a one-dimensional dimerized optical superlattice, thus 
obtaining 

#bdi = 2(i - 5)4„-iC 2 » + (t + 5)c\ n c 2n -i + H.c, (17) 

n 

where 6 quantifies the different tunneling strength between superlattice sites |35j . On 
the other hand, the Creutz ladder is described by 

H Am = Ke- ie al +1 a n + Ke l6 b\ l+1 b n + Kb\ +1 a n 

n 

+ Ka\ l+1 b n + Ma) n b n + H.c, (18) 

where K, M are tunneling strengths, and 8 is a magnetic flux piercing the ladder. 
This requires two Zeeman sublevels to be assigned to the fermion species a n ,b n , and 
a one-dimensional laser- assisted tunneling U ai = diag{e , e }, U &1 — ia x , together 
with the Raman on-site operator of strength M 63 . 

It is possible to continue this approach, proceeding thus to higher dimensions 
and different topological classes. Prominent examples would be the honeycomb time- 
reversal breaking Haldane model [Ml ES| for the D — 2 topological insulator in class 
A, or the time-reversal Kane-Mele model in the honeycomb lattice in class All [6"6"] . or 
other optical-lattice geometries [27j EH [68] • Rather than following this route, we shall 
explore a different approach that is better suited to the superlattice-based simulator 
introduced above. Indeed, we shall argue that this quantum simulator allows the 
reproduction of most the topological phases in table [5[ 



Dimensional- Reduction Approach. In this case, the starting point is the quantum 
simulator of Z?-dimensional Kaplan fermions in equation (16). Depending on the 
particular choice of Dirac matrices, the inverted-mass regime shall correspond to a 
different class of topological insulators. Besides, in some situations, a dimensional 
reduction [53 that amounts to the increase of the optical-lattice depth in one direction, 
connect us to a different lower-dimensional class. We rewrite the full Hamiltonian 



= J2*Hp<i)H£Mpd), #Di(p d ) = 

d.Pd 



Pd + TOK, 



(19) 



to — ^2 v {—^) d "rnv, and where the Dirac matrices 
ce~,p shall be selected so that the T,C,S symmetries are explicitly broken or 

Hamiltonian preserves these symmetries 



where the Wilson mass is 

„d 



preserved. This translationally-invariant 
when the following conditions are met 

T: 4 [ffgjf-pd)]' 



C : 
S 



9 



(-Pd)r 



6 T = +ffgi(Pd) 

e c = 



^m(Pd); 



e 



ekffg I ( Pd )e c e; 



(Pd); 



(20) 
(21) 

(22) 
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Table 6. Quantum simulator of Topological Insulat ors, We list different 
realizations of Wilson-fermion Hamiltonians in equation ^19\ that directly lead 
to several classes of topological insulators. Each class, characterized by the 
discrete symmetries T,C,S, where T 2 = ©t©t = -t^' anc ^ ^ = ®C©c = 
Besides, each class has a Wilson-fermion representative with a particular choice 
of the Clifford algebra «„,/3 that depends on the dimensionality D and the 
corresponding symmetries. We also highlight the topological insulators that 
can be obtained by dimensional reduction from a parent Hamiltonian, such as 
All, D = 3 =-> All, D = 2, or A, D = 2 ^ AIII, D = 1. We also list the unitary 
matrices 9t, @c involved in the definition of the discrete symmetries. 



Class 


D 




OLy 


OL z 


p 


T 


©c 


r 


c 


s 


CII 


3 


a z <g o x 


o z <g Oy 


O z <g o z 


o z <g I 2 


il (g) CTy 


io x (g (7 j, 


-i 


-i 


1 


AIII 


3 


o z <g o x 


Z <g Oy 


o z (g o z 


Oy (g I 2 


il (g (Ty 


ifjj; (g (Tj, 








1 


All 


3 


o z <g> o~ x 


Z <g Oy 


O z <g o z 


cr^ (g)I 2 


il (g CTy 


idj; (g (Tj, 


-i 








^ All 


2 


o z <g o~ x 


Z ® Oy 




CTa; (gl 2 


il (g CTy 


io x (g (Tj, 


-i 








A 


2 


0~x 


Oy 






I 


i^x 











AIII 


1 


0~x 






<7z 


I 


i^x 








1 


CII 


1 


a z <g o x 






o z <g I 2 


H <g CT y 


\O x (g (Tj, 


-i 


-1 


1 



where Ot, ©c are some unitary matrices. In table[6j we list the symmetry properties of 
different Kaplan-fermion Hamiltonians. It is important to note that these symmetries 
might correspond to the exact symmetries in nature (e.g. when considering the 
hyperfme levels {\F, mp), \F, — mp )}, the time-reversal symmetry given by 6t = 
io y exactly correspond to time-reversal symmetry in nature). Conversely, these 
symmetries might otherwise correspond to the algebraic properties of the effective 
Hamiltonian. Let us emphasize, however, that as far as the disorder respects such 
symmetries, the robustness of the edge excitations is guaranteed. It would be of 
the greatest interest to design disorder breaking or preserving such symmetries, 
generalizing the studies on Anderson localization [7] . 

In table [6j we have listed the different topological insulators that can be simulated 
with our scheme. As shown in [29] for the particular case of three-dimensional AH 
insulators, the laser parameters can be controlled so that an odd number of Wilson 
masses are inverted. This mass-inversion occurs through a gap-closing point, and 
thus a quantum phase transition between a normal band insulator, and a topological 
one occurs. This new phase is characterized by an odd number of massless fermionic 
excitations (i.e. massless Dirac fermions) bound to the boundaries of the system, and 
protected by a topological invariant. In the three-dimensional case, this corresponds 
to an axion term that modifies the response of the system according to the so-called 
axion electrodynamics |69j . Remarkably, table [6] contains all the relevant information 
to explore the exotic properties of different topological insulators in a superlattice 
based experiment with ultracold atoms. 
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6. Conclusions 

In this article, we have presented a concrete proposal for the realization of laser- 
assisted tunneling in a spin-independent optical lattice trapping a multi-spin atomic 
gas. Remarkably enough, it is possible to tailor a wide range of spin-flipping hopping 
operators, which opens an interesting route to push the experiments beyond the 
standard superfluid - Mott insulator transition. The scheme we have presented 
combines bi-chromatic lattices and Raman transfers, to adiabatically eliminate the 
states trapped in the middle of each lattice link. These states act as simple spectators 
that allow us to assist the tunneling of atoms between the main minima of the optical 
lattice. This mechanism is clearly supported by our numerical simulations of the time- 
evolution of the atomic population between the different optical-lattice sites. Even if 
we focus on fermionic 40 K, we stress that the ingredients of this proposal do not rely 
on the atomic statistics, and could be thus used for all the alkalis. 

We believe that such a device could have important applications in the quantum 
simulation of non-interacting lattice field theories, which are characterized, in their 
discrete version, by on-site and nearest-neighbor hopping Hamiltonians. Once the 
fields of the theory to be simulated are mapped into the atomic hyperfine states, the 
desired operators correspond to population transfers between such levels. The former 
can be realized by standard microwaves, whereas the latter might be tailored with the 
laser-assisted schemes described here. 

Even though interactions are at the heart of a plethora of interesting effects, 
non-interacting fermionic theories already encompass a number of phenomena whose 
experimental realization would be of the greatest interest. In the second part of the 
article, we analyzed interesting physical models which could benefit from our proposal. 
In particular, we focused on relativistic field theories, and showed that there is a zoo 
of relativistic lattice fermions that can be addressed with this platform. Besides, we 
presented a toolbox to design particular assisted tunneling processes that lead us to 
the physics of topological insulators. Remarkably enough, this quantum simulator 
turns out to be extremely versatile, since most of the phases of the periodic table of 
topological insulators can be addressed. 

Finally, let us comment on the possible combination of this proposal with the 
control of interactions already achieved in cold-atom gases. This might eventually 
boost experiments into regimes where classical numerical simulations fail, which we 
leave here as an outlook for future work. In particular, the problem of robustness of 
topological orders (classified for non-interacting theories) with respect to interactions 
is one of the most important challenges of the modern condensed matter TOJ . We 
believe that a direct combination of our setup with Feshbach resonances will provide 
important insights into this unsolved question. 
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